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1 INTRODUCTION 



According to present experimental constraints, the universe, nowadays, is almost spatially 
flat; moreover, the effect of the radiation term, in the Friedman's equation, is also neglige- 
ble. However, one should distinguish between the description of the universe, as it appears 
"now" ,and the description the description of how it was and how it will be. 

Most papers dealing with the determination of the cosmological parameters, nowadays, 
deal with the universe as it appears "experimentally" , now. Many other articles deal with 
the complicated issues related with the mechanism(s) of inflation, reheating etc. 

In the present article, starting with recent experimental constraints (|6|,0,Q,|^]) on 
the density parameters ^Im, ^A, we want to study analytically the whole history of the 
universe, i.e., how it was {after the first few minutes) and how it wil be. A detailed analytic 
discussion of Friedman universes with cosmological constant and radiation pressure was 
performed years ago j^), but at that time, the cosmological constant was believed 
to be zero by many, so that the above work - although justified by a particular phe- 
nomenological analysis ( |jl^) - could be considered only as only a purely mathematical 
exercise. 

At the light of recent experimental results, part of the work carried out in (Q, [^) 
becomes phenomenologically relevant : we do not plan to give a thorough discussion of 
analytical solutions in all possible cases, since this was done already by us, but restrict 
our attention to those cases that are compatible with recent observations and present an 
updated analytic discussion. 

We are not concerned, in this paper, with the problem of inflation (there are several 
theoretical scenarios ... ) as it is clear that, at such early periods, Friedman's equation 
is anyway not expected to provide a good mathematical model of the evolution of the 
universe. It looks however reasonable to suppose that there was a time, in some remote 
past, where the universe was already almost homogeneous and still rather hot. One can 
then use Friedman's equation to extrapolate, back in time, to this early epoch. Putting 
artificially the radiation term to zero is an unnecessary simplification that would prevent 
one to perform useful extrapolation towards this remote past. One can actually keep 
that term at no cost since cosmological models with a non-zero cosmological constant and 
radiation pressure can be studied analytically. Moreover, the experimental constraints 
on the radiation term are quite well known thanks to the measurement of the Cosmic 
Microwave Background radiation (CMB). 

A similar comment can be made about the value of the constant parameter (k = ±1 or 
0) that allows one to distinguish between spatially closed {k = 1) and open (A; = or — 1) 
universes: it may be that the reduced cosmological curvature density $7^ is very small, 
nowadays, and experimentally compatible with zero; however, this quantity is function of 
time, and setting artifically the constant A; to prevents one from studying some possibly 
interesting features in the history of the universe (for instance, the existence of an inflexion 
point at some particular time during the expansion, requires k = 1). 

For the above reasons we shall not drop, in general, the radiation term and will keep 
k as an unknown discrete parameter with values ±1 or 0. Our main purpose, however, 
is to study the effects of a non zero (and positive) cosmological constant on the analytic 
solutions describing the dynamics of our universe and on the large scale geometry. 

In the coming section, we shall recall the relations between different cosmological 
quantities of interest. One should remark that experimentally available parameters do not 
necessary coincide with the quantities that provide a good mathematical description for 
the solutions of Friedman's equation. 
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In the following section we shall start with this equation (written with a "good" set 
of variables), and shall analyse possible histories of our universe, both qualitatively (using 
a well known associated mechanical model |jl^) and analytically (using mathematical 
techniques (that can be traced back to Q) explained in Q and used in Q, [^). 

Finally, in the last section, we investigate the influence of the selection of a cosmo- 
logical model matters, on "changes of redshift charts". We answer the following question 
(following again Q): if some observer, located somewhere in the universe, plots the red- 
shift of every celestial object against its direction, what is the redshift that would be 
measured, for the same objects, by another observer, at the same moment of time (we 
have an homogeneous cosmology!) but situated somewhere else - and possibly very far - 
in the spatial universe ? These transformation laws involve not only the geometry of our 
universe (relative location of the two points) but also its dynamics (in particular the value 
of the cosmological constant). 

If the distribution of matter in the universe posesses some particular large scale features 
(for instance a "symmetry" of some kind), our solution of the above problem becomes 
phenomenologically relevant. Indeed, there is no reason to believe that our position (our 
galaxy) is ideally located with respect to such putative elements of symmetry, and it 
would therefore be necessary to perform a change of position to recognize such features. 
For instance, the absence of matter in a three dimensional shell centered around some 
particular point would not appear as such when analysed from a galaxy situated far away 
from this special point (this was actually the hypothesis made in [p!o|]). 

This paper presents two very different results: one is concerned with analytical solu- 
tions describing the dynamics of the universe, the other is concerned with large scale (and 
non-euclidean) geometry of our 3-dimensional space manifold. In both cases we want to 
draw the attention of the reader on the novel features presented here. 

When the cosmological constant A is zero, the dynamics of our universe can be de- 
scribed in terms of elementary trigonometric (circular or hyperbolic) functions. When A 
is non zero, the corresponding Friedmann's differential equations, cannot be expressed in 
terms of trigonometric functions. The remark we make in this paper, (and in an earlier 
one (Q)) is that that this equation may be handled with a somewhat more general class 
of periodic functions. 

Different explicit expressions can be found in papers by G. Lemaitre ([l[), but these 
expressions, involving elliptic integrals, are not very transparent when written as a function 
of cosmic time. 

However, a simple trick -namely, writing Friedmann's differential equation for the 
inverse "radius of the universe " in terms of conformal time rather than for the radius 
in terms of cosmic time- allows one to find explicit, analytical solutions, in terms of very 
basic elliptic functions, namely, the Weierstrass elliptic function that have been arounds 
for well over a century. We have not been able to find any description of this simple trick 
in various books or review articles (e.g. [^2|) - even the most recent ones like for instance 
(lH]) " devoted to the description of a universe with cosmological constant. 

We should maybe stress the fact that the Weierstrass elliptic functions are very simple 
(if not elementary!) objects, that generalize both circular and hyperbolic trigonometric 
functions; for example, the Weierstrass V function is a doubly periodic function in the 
complex plane (in the limit when one of the periods becomes infinite, it can be expressed 
in terms of usual trigonometric functions). It also satisfies interesting relations, such as a 
duplication formula analogous to the one valid for trigonometric functions. This leads to 
a relation between cosmological quantities measured at different times (more about this 
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in section 4.2.2| ). 

Another subject, mostly disconnected from the problem of finding explicit analytical 
solutions, is the description of large scale geometry of our three-dimensional (spatial) uni- 
verse; this aspect of our work was already presented in the above introduction. However, 
in order to dissipate possible misunderstandings, let us stress again the fact that the cor- 
responding geometry involves the curvature of a 3-dimensional curved manifold (not like 
the 2-dimensional sphere which is the object of spherical trigonometry), so that even if our 
formulae (concerning e.g. arc-length of the sides of a geodesic triangle) are reminiscent 
of usual spherical trigonometry, they are actually quite different. This part of our work 
can be considered as a physically relevant but not totally trivial mathematical exercise 
in riemannian geometry. As a matter of fact, the corresponding formulae do not seem to 
be available in mathematical treaties of differential geometry (probably because they are 
"nothing else" than a particular example of calculations that one can do on any riemannian 
manifold when its underlying space has the topology of a Lie group). 

From the (astro) physical point of view, we do not believe that such formulae have 
been presented elsewhere. 
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2 RELATIONS BETWEEN COSMOLOGICAL PARAM- 
ETERS 



2.1 Geometrical Background for Friedman-Lemaitre Models 

Most of the present section is well known and can be found in standard textbooks. The 
dimensionless temperature T(t) which seems preferable to the radius R{t)) if one wants to 
perform an analytic discussion of Friedman's equation, was introduced by Q. Our purpose 
here is mostly to provide a summary of useful formulae and present our notations. 

2.1.1 Units 

We chall work in a system of units in which both h and c are set to 1. All quantities with 
dimention are thus given as powers of length [cm]. 

In our units, energies are in [cm~^]. Since temperatures (to be denoted by T in this 
paper) appear only in products ksT (where fcs is the Boltzmann constant), we can take 
take the temperature as dimensionless, and ks as an energy. 

2.1.2 Metric and Topological Considerations 

As it is well known, these models assume that the universe is, in first approximation, 
homogeneous and isotropic. 

The metric There are three possibilities: 

• k = l,ds^ = -dt^ + R^[dx^ + sin^ xidO"^ + sin^ 9d(t?)] 

• k = -l, ds^ = -dt^ + R^[dx^ + sinh^ xidO"^ + sm^ Odcf^)] 

• A: = 0, (is2 = -dt^ + R^[dx^ + x^idO"^ + sin^ Odcf^)] 

The topology The first case /c = 1, is sometimes called the "compact case" (closed 
universe), since the spatial universe is then compact; the global topology for space-time 
can be taken as S*^ x R but one should remember that the above expression for the metric 
is only local, so that one could take, in place of S^, any quotient S^/P of the 3-sphere by 
a finite group F operating without fixed point (for instance one of the binary tetrahedral, 
octahedral or icosahedral subgroups of SU (2)). In the case of 5^, the full isometry group is 
SU{2) X SU{2) since itself is homeomorphic with the group of unimodular and unitary 
matrices. 

In cases k = Q oi k = —1, the space component is non compact (open) and the universe 
is respectively called "flat" or "hyperbolic". 

On intuitive grounds, and independently of experimental fits, one could argue that the 
model with A: = 1 is preferable to the others since Friedman's equation was derived by 
assuming homogeneity of the stress-energy, something that is hard to achieve in an open 
universe if the total quantity of matter is finite. 

2.1.3 Friedman's Equation 

This equation is obtained by writing Einstein's equation G^^ = SvrT^jy where T^y is the 
stress energy tensor of a perfect fluid and where G^y is the Einstein tensor associated with 
the above metrics. 

The average local energy splits into three parts : 
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• Vacuum contribution p^ac = where A is the cosmological constant, 

• Radiation contribution prad such that R'^{t)prad{t) = const. = g^Cr, 

• Averaged matter contribution pm such that R^{t)pm{t) = const. = Q§QCm , 
The evolution of R{t) is governed by the Friedman's equation 

1 , dR , r, Cy C'lri k A 

— f — r = — H — 1 

R'^^ dt^ R^ R3 R'^ 3 

It is convenient to introduce the conformal time r, defined, up to an additive constant by 

dt 

'''=R 

This variable is natural, both geometrically (it gives three-dimensional geodesic distances) 
and analytically (as we shall see later). The quantities C^, Cm and A have dimensions LP', 
L, L^^ respectively. In order to replace i? by a dimensionless quantity, we introduce the 
characteristic length scale of matter 

A 

Ar 



9C2 



and replace R{t) by the reduced dimensionless temperature 



kI"r{t 

With this new variable, Friedman's equation becomes 



= . . (1) 



{dT/drf = aT^ + -T^ - kT^ + - (2) 
3 3 

where 

Ac 

and 

a = Cr^c 

are two constant dimensionless parameters. 

The reason for calling T the "reduced dimensionless temperature" is that it is propor- 
tionnal to the temperature T of the black body radiation. Indeed, prad = ^cfT^ where a 
is the Stefan-Boltzmann constant. Since prad = o-^aAcT"^, one finds 



SttG 

3 aA, 



SvrG 4c7 



(3) 



This may be the right place to recall that Lemaitre (||T|) did, long ago, an analytic 
study of the solutions of Friedman's equations; his discussion, made in terms of R and 
t involves, of course, elliptic functions, but it is only when we express the cosmological 
quantities in terms of T{t) (the reduced temperature, as a function of conformal time) that 
these quantities can be written themselves as elliptic functions with respect to a particular 
lattice (the associated Weierstrass V function appears then naturally) . The fact that such 
a direct link with the theory of Weierstrass can be established should be clear from the 
fact that, with our parametrization, the right hand side of equation § is a polynomial of 
degree four. 
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The Hubble function decribing the rate of expansion is defined, as usual, by H 
R-^{dR/dt) and can be written, in terms of T(r) as 



1 dT 



Another useful quantity is the deceleration function q = —R{dR/dt) '^{d'^R/dt^) that 
can be expressed in terms of T{t) as follows: 



Actually, the choice of the sign (and of the name) turns out to be a historical mistake, 
since the recent bounds on the cosmological constant lead to a negative value for q (an 
accelerating expanding universe). 

By multiplying Friedman's equation by Ac/H"^-, one obtains the famous relation : 



l = VLr + Vlm + ^k + ^lv (4) 

with 

= -kT^^ (5) 

^r. = ^T^Al (6) 

= aT^^ (8) 

Notice that the famous "cosmic triangle" of - a kind of cosmic Dalitz' plot, becomes 
a triangle only if one decides to forget the radiation term (then, indeed, 1 = Qm+^k + ^A)- 
In terms of the usual temperature T, the radiation terms reads: 

2.1.4 Description of the History of the Universe by the Reduced Temperature 
Function T{t) 

• Rather than using R{t), we describe the history of the universe by the function 
T{t). This function can be intuitively thought of as a dimensionless inverse radius; it 
is proportionnal to the usual temperature T. Because of the phenomenon of expan- 
sion, these quantities T and T decrease with r. The argument r itself, a dimensionless 
arc length, measures time : for instance, if the spatial universe is (hyper)spherical, 
r is a measure (in radians) of the length of the arc spanned by a photon that was 
produced at the big bang. A given solution (a given "history" of a - dimensionless 
- universe) is fully characterized by the solution of a differential equation depending 
on the two dimensionless parameters a and A. Of course, one has also to specify 
some initial value data (one takes T(r) — > oo when r ^ 0). As we shall see in a later 
section, T(r) is a particular elliptic function; the history of our universe is then also 
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fully encoded by the two Weierstrass invariants (§2,93), or, equivalently, by the two 
periods loi,uj2 of this elliptic function. The Weierstrass invariants are given, in terms 
of parameters k, a and A, by the algebraic relations 



= 12 + T- ^''^ 



1 , n , , aXk 
9s = -,ik'-2X)-^. (11) 

The formulae giving the two periods uji , 102 in terms of the cosmological parameters 
A;, a and A (or 31,52) involve elliptic integrals and will be given later. 

• Dimensional quantities are obtained from T(r) thanks to the measurement of a single 
function having dimensions of a length. Such a function is usually the Hubble function 
H = H{t). The behaviour of R{t) is given by the two parametric equations R{t) and 
t(r). As we shall see later, in all cases of physical interest, t(r) reaches a logarithmic 
singularity for a finite conformal time t/ (the universe expands forever, but as t — 00 
the arc length associated with the path of a single photon goes to a finite value rj). 

• Finally, it remains to know "when"' we are, i.e., the age of the universe. This can 
be expressed in terms of the dimensionless quantity To (a particular value of r) or, 
more conventionally, in terms of to = t{To). Experimentally, one measures the Hubble 
"constant" Ho = H{to), i.e., the value of the function H{t), now. 

2.1.5 Cosmological Quantities 

We just give here a list of most cosmological quantities of interest. 

Dimensionless Quantities We have the constant parameters a, X, k = 0,±1, the 
time-dependent dimensionless densities r2fc(r), J1a(t), f^m (''")) ^riT) and the reduced tem- 
perature function T{t). Of interest also is the deceleration function g(r). Of course, the 
conformal time r itself is a non-constant (!) dimensionless quantity. The black-body tem- 
perature T (temperature of the blackbody microwave radiation, experimentally expressed 
in degree Kelvin ) is itself dimensionless - sec our remark at the beginning of this section. 
It is proportional to the dimensionless reduced temperature T (or to the inverse radius 
R). 

Dimensional quantities 

• The Hubble function H{t). This quantity is usually chosen to fix the length scale. 

• The critical length scale A^. It can be thought of as giving a measure of the "total 
mass" of the universe. More precisely, if the universe is spatially closed and has 
the topology and metric of a 3-sphere, its total mass is M = Aii^R^pm so that 

■^c = (2^) • 

• The radius R. Intuitivelly, it is a measure of the "mesh" of our spatial universe. 

• The age t of the universe (cosmic time) . 
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2.1.6 Expression of the Cosmological Quantities in Terms of the Hubble Func- 
tion H and of the Dimensionless Densities J^a, fifc, J^m, Jl^ 

Since most experimental results are expressed in terms of the dimensionless densities 
J7a, fifc, 0^, we express all other cosmological quantities of interest in terms of them. 
Typically, dimensional quantities also involve the value of the Hubble function H{t). All 
these formulae can be obtained by straightforward algebraic manipulations. 

Constant dimensionless quantities (parameters) 



Time-dependent dimensionless quantities (assuming k ^G) 



2 2 \Q,k\ 

q=^ -nK + 0.r = 3^ +nk + 2ar-\ (15) 
The temperature of the CMB radiation: 

87rG4(7 

The conformal time r itself can be found numerically, once a value of T is known, by 
solving the (non algebraic) equation expressing T as a function of r (see the next section). 

Constant dimensional quantities (parameters) 

A = ^H'^^k (16) 



A. = 1M^^ (17) 

^ ^777, 



Time-dependent dimensional quantities (assuming k ^ Besides the Hubble 
"constant" H itself, one has to consider also the radius 

TAc^/^ 

the cosmic time (the age of the universe corresponding to conformal time r) given by 
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2.2 Dimensions of the Cosmological Quantities 

In our system of units, all quantities are either dimensionless or have a dimension which 
is some power of a lengh ([cm]). We gather all the relevant information as follows: 

t rsj R ~ [cm] 

Gr^Cr^ [cm^] 
H ^ ks ^ energy ~ [cm~^] 
A ~ A. ~ fcm"^l 



-4i 



Pvac ~ Prad ^ pm ^ ^ [cm 

Here a denotes the Stefan constant (remember that degrees Kelvin are dimensionless). 
Finally we list the dimensionless quantities: 

~ i^A ~ ^ Q ^ T T ^ T ^ a ^ X ^ [cmP = 1] 
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3 EXPERIMENTAL CONSTRAINTS 



3.1 Experimental Constraints from High-redshift Supernovae and Cos- 
mic Microwave Background Anisotropics 

In general we shall add an index o to refer to present values of time-dependent quanti- 
ties (like Ho, ni, n"r, ^m, ^1, 

Qo, To or To). Remember that k,a,X,A,Ac are constant 

parameters. 

The analysis given by HI (based on experimental results on the latest cosmic microwave 
background anisotropy and on the distant Type la supernovae data |^ gives 0.13 < 
Q'^ < 0.43 and +0A0 < 0^ < +0.80. In order to stay on the safe side, we do not 
restrict ourselves to the more severe bounds that one can obtain by performing a combined 
likelihood analysis that would use both sets of data coming from anisotropy of CMB 
radiation and distant type la supernovae. 

Moreover, many experiments, nowadays, give a present value of the Hubble function 
close to Ho = h 100 km sec~^Mpc~^, with h ~ 0.66, so that H~^ ~ 14.2 x 10^ yr = 
13.25 X 10^7 cm. This is the value which we shall use in our estimations. 



3.2 Implications of Experimental Constraints on Friedman-Lemaitre Pa- 
rameters 

With h = c = 1, the value of the Stefan-Boltzmann constant is a = — gjp- = 59.8cm , ks 
is the Boltzmann constant and the Newton constant is G = 2.61 x 10~^^ cm^. Using the 
above value for Ho and formula |9| (with To = 2.73K), we obtain Q° ~ 5.10 x 10~^. 

Using now the "cosmic triangle relation" (equation ^) together with the experimental 
bounds on Q'^ and previously recalled, one finds —0.23 < < 0.47. Remember that, 
with the present conventions, 0^ < when k = +1. 

The formula 16 given in the last section leads to +0.68 x 10~^^ cm^^ < A < +1.36 x 
10-56 cm-2. 

In order to obtain good estimates for the remaining quantities, one would need a 
better measurement of the curvature density fi^. Indeed, formulae |l^, |l^, together with 
the previously given bounds on 0^ only imply the following for the dimensionless constant 
parameters: a < 6.3 x 10^^ and A > 0.44. 

Notice that A is quite sensitive to the independent values of 17$^ and (indeed, 
fifc, appearing in the denominator of equation 12, can be very small); for instance, the 
values = 0.13, ni = 0.40 lead to A = 0.44 but = 0.13, = 0.80 lead to A = 
266.68. By the way, one should stress the fact that the behaviour of analytic solutions of 
Friedman's equations is essentially determined by the value of the (dimensionless) reduced 
cosmological constant A, which can be rather "big", even when the genuine cosmological 
constant A (which is a dimensional quantity) is itself very "small" . 

As already mentionned, present experimental constraints only imply A > 0.44, but 
one should remember that ~ 1 — — should be negative for a spatially closed 
universe {k = +1); therefore, if, on top of experimental constraints on Qrn and we 
assume that we live in a spatially closed universe (a -reasonable- hypothesis that we shall 
make later, for illustration purposes) then 0^ + ^^A > 1 and the constraint on A, as given 
by equation 12 is more tight; an elementary variational calculation shows then that the 
smallest possible value of A, taking into account both the experimental constraints and 
the hypothesis k = +1, is obtained when both flm and Oa saturate their experimental 
bounds (values respectively equal to 0.43 and 0.8), so that fi^ > —0.23 and A > 82. 
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One should not think that the value of A could be arbitrary large: we shall see a little 
later (next subsection) that, for experimental reasons, it has also to be bounded from 
above (condition A_ < A < A+). 

Using equations 14, one obtains easily the bounds To > 0.41, —0.75 < qo < —0.18 



and Ac < 1.6 x 10~^^cm~^ for the present values of the reduced temperature T and 
deceleration function q and for the critical length scale Ac. 

In view of the experimental results (17° small and fi^ compatible with zero), one maybe 
tempted to make the simplifying hypothesis 1 = OJ^ + J^^, i.e., set both terms il^ and to 
zero in the relation 1 = $7^, + ^k + + ^r] notice however that assuming the vanishing of 
ilfc + fir at all times (without assuming the vanishing of each term, individually) is totally 
impossible, as these densities are not constant and it is easy to see, from the definition of 
these quantities, that such a relation can only hold at at a single moment; it is obvious 
that the radiation contribution, 0°, although small, is not strictly zero. Moreover, putting 
artificially the constant k to zero is certainly compatible with the present experimental 
data, but one should be aware of the fact that the curvature density fl^ is not a constant 
quantity and that setting it to zero at all times is an artificial simplification that was 
probably not justified when the universe was younger . . . 

3.3 A Curious Coincidence 

The parameters Cm and Cr appearing in Friedman's equation, and measuring respectively 
the matter contribution and radiation contribution to the average local energy, are, a 
priori , independent quantities. For this reason, the relation between T (the temperature 
of the cosmological background radiation) and T (the reduced temperature) involves Cm 
and Cr- 



r 24 

aAc^4 _ 3 '-"'•Fc 



SvrG Aa SvrG Aa 

One may consider he special case of a universe for which the two functions T(r) and 
T(r) coincide; this amounts setting 

Ct 60 10 ^ 

In this case, the present values of the densities fi^ and 0^ are no longer independent: 
r2° = 5.1 X 10~^ is, as usual, exactly known since To = 2.73K is known; again, as usual, 
we have flm = fl^fcl^'o but if we measure (for instance) and assume To = To, we find 
the value of Choosing then, for instance, /c = +1 (a closed spatial universe), so that 
ni = -\ni\, we deduce from the equation ^1 + + + fl'^ = 1. 

We do not see any theoretical reason that would justify to set T = T, since these two 
quantities could very well be different by several orders of magnitude (even with the rather 
small value of the radiation density 17°), nevertheless ... by a curious coincidence, if one 
takes the two densities and 17^ to be numerically equal to their highest possible values 
compatible with present experimental bounds (0.43 and 0.8), one finds a value of To equal 
to 2.80; this is very near the experimental value 2.73 of the temperature To- We have no 
explanation for this "cosmological miracle". 
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4 ANALYTIC BEHAVIOUR OF SOLUTIONS 



4.1 Qualitative Behaviour of Solutions 

Equation (|2|) can be written 

with 

Vo^,k{T) = -aT^-^T^ + kT^ 

This is the equation of a one-dimensional mechanical system with "coordinate" T, potential 
V^^fc(T) (shown in Fig |^) and total energy A/3. The kinetic energy beeing non negative, 
the associated mechanical system describes a horizontal line in the (y{T),T) plane but 
never penetrates under the curve V^,fc(r) (that would correspond to r imaginary). The 
length of the vertical line segment between a point belonging to the curve and a point 
with same value of T but belonging to the horizontal line A/3 (described by the associated 
mechanical system) is a measure of (^)^- 

For a given value of a (the radiation parameter), and for A: = ±1, the curve Va,k{T) has 
typically two bumps (two local maxima). For k = 1 (closed universe), the right maximum 
occurs for a positive value of T and V{T) whereas, for k = —1 (hyperbolic case), this 
maximum is shifted to T = and V{T) = 0. For k = (flat case), the right maximum 
disappears and we are left with an inflection point at T = 0, V{T) = 0. 

In this paper, we almost always suppose that a is non zero; however, one should notice 
that if a = 0, the curve V{T) becomes a cubic (Fig. |2|) and the left hand side maximum 
disappears (it moves to —00 as a goes to 0). 

Notice that figure || describes only the qualitative features of our space-time history, 
since, for reasonable values of a and A (i.e. values compatible with experimental con- 
straints), the vertical coordinate of the left hand side maximum should be at least 1000 
times higher than the vertical coordinate of the right hand side maximum. 

The experimental constraints also show not only that A is positive (since is) but that 
it is bigger than A_/3 (the vertical coordinate of the right maximum) and, at the same time 
much smaller than A+/3 (the vertical coordinate of the left maximum). Indeed, for small 
values of a, the positions of the extrema (cf. next section) are given by A_ ~ k(l — 3ak) 
and A+ ~ i6^(l + 12A:a); the previous bounds on A show then directly that A > A_; 
moreover, the inequality A < A_|_ is equivalent to < and this clearly holds. 

Unfortunately, this does not lead to a stringent constraint on A istself since the bound 
Q < 6.3 X 10~^ only implies A+ > 2.5 x 10^. 

The associated mechanical system is therefore such that it describes an horizontal line 
like the one displayed on Fig. |l[ Typically, a given universe starts from the right of the 
picture (T — > 00 corresponding to the big bang), moves to the left (slowing down), until 
it reaches the vertical axis (T = 0). This actually takes place in a finite conformal time 
Tf but, as we shall see, it corresponds to a cosmic time t 00, so that, for the universe 
in which we live, history stops there. However, the solution can be continued for T < (a 
negative radius R) until the system bumps against V{T) to come back to right infinity; 
the system then jumps to left infinity, follows the horizontal line (to the right) till it bumps 
on V{T) again, and comes back. This round trip of the associated mechanical system is 
done in a (conformal) time 20;^ ~ a period of the corresponding elliptic function. Let us 
stress the fact that only the first part of the motion (from right infinity to the intersection 
with the vertical axis) is physically relevant for ordinary matter. 
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Figure 1: Potential for the associated mechanical system. Case a 7^ 
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Figure 2: Potential for the associated mechanical system. Case a = 

Since the radius R is proportionnal to 1/T (see equation ||), the discussion, in terms 
of R is of course different: the system starts with R = (big bang), and expands forever; 
moreover the expansion speeds up in all three cases k = ±1, A; = 0. The only particularity 
of the case A; = 1 is that there is an inflexion point (coming from the existence of a positive 
right hand side maximum for the curve V(T)) : the expansion speeds up anyway, but there 
is a time tj for which the rate of expansion vanishes. 

4.2 The Elliptic Curve Associated with a Given Cosmology 
4.2.1 General Features 

Let us call Q{T) the fourth degree polynomial that appears at the right hand side of 
Friedman's equation (eq. |2|). Let Tj be any one of the (possibly complex) roots of the 
equation Q{T) = 0, then the fractional linear transformation 



4 T-Tj 24 



A'^=P(y)=4y^-g^y-g, (21) 



brings to the form 



where 52 and are the two Weierstrass invariants given in eqs. ( [IC|JTl| ). The solution to 
the previous equation is well known: y = V{t; §2, gs) where V is the elliptic Weierstrass 
function associated with the invariants 52 and (73. One should remember that, given a 
lattice L in the complex plane, an elliptic function / with respect to L is non constant 
meromorphic function that is bi-periodic, with respect to L (so that /(t + u) = /(t) 
whenever u belongs to the lattice L). It is known (since Liouville) that if a is an arbitrary 
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Figure 3: An elementary periodicity cell 



complex number (including infinity), the number of solutions of the equation /(r) = a is 
independent of a, if multiplicities are properly counted; this number is called the order 
of f. It is useful to know that any rational function of an elliptic function is also elliptic 
(with respect to the same lattice) but that its order will coincide with the order of / only 
if the transformation is fractional linear (like the transformation (|20|)). Finally, one should 
remember that the order of an elliptic functions is at least two and that the Weierstrass V 
function corresponding to a given lattice is defined as the elliptic function of order 2 that 
has a pole of order 2 at the origin (and consequently at all other points of L) and is such 
that — V{t) vanishes at r = 0. The two elementary periods generating the lattice L 
can be expressed in terms of §2 and gs, but conversely, the Weierstrass invariants can be 
expressed in terms of the elementary periods of the lattice L. 

These old theorems of analysis lead directly to the fact that T{t) is an elliptic function 
of order 2 and that other quantities that are rational functions of T (like the Hubble 
function H, the deceleration function q or all the cosmological densities ^Im, 0^,0^, J^a) 
are also elliptic with respect to the same lattice (but they are not of order 2). In this sense, 
one can say that our universe is fully described by the elliptic curve = — g2y — gs 
characterized by the two Weierstrass invariants g2, gs, or, alternatively, by two elementary 
periods generating the lattice L. 

In this (physical) case, the elementary periodicity cell is a rhombus (figure ^) with 
vertices {0, Ur — iiOi, 2L0r,^r + iuJi}. Another standard notations for periods are 2uj2 = 2uJr, 
—2uJ3 = LOr + i^i and 

— 2uii = LOr — iuJi 

(so that (jji + i02 + UJ3 = 0). The transformation expressing the half-periods a'1,2,3 in terms 
of the Weierstrass invariants 52 and 53 (or conversely) can be obtained, either from a direct 
numerical evaluation of elliptic integrals (see below) or from fast algorithms described in 



I nil; one can also use the Mathematica function WeierstrassHalfPeriods[{g2,g3}] as 
well as the opposite transformation Weierstrasslnvariants[{ivi,uj2}]- In the case a = 0, 
the value t = ujr turns out to be a global minimum of T{t), therefore another possibility 
to get uJr in that case is to find numerically the first zero of the equation dT{T)/dT = 0. 



4.2.2 Neglecting Radiation : the Case a = 

When a = 0, the two poles of T(r) coincide, so that the dimensionless temperature, up 
to shift and rescaling, coincides with the Weierstrass "P-function iself (the fractional linear 
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transformation 2C takes the simple form T = 6y + k/2). Notice that V is an even function 
of r. 

In that case 

T{T) = 6[V{T;g2,g3) + k/12] (22) 

A typical plot of T(r) for experimentally reasonable values of the cosmological parameters 
is given in figures ^ (closed universe: A; = 1) and |5| (open universe: A; = — 1). 

In both cases, the physically relevant part of the curve is given by the interval < r < 
Tf since r = corresponds to the big bang (infinite temperature) and r = t/ to the end of 
conformal time: the cosmic time t = i(r) reaches a logarithmic singularity when r Tf, 
so that when t — > oo, the arc length described by a photon born with the big bang tends 
toward the finite value Tf, for values of A compatible with recent observations, this limit 
is strictly less than 2tt so that it will never be possible to see the "back of our head" , even 
if the universe is closed and if we wait infinitely many years .... 

Notice that both curves show the existence of an inflexion point denoted by tj on the 
graphs; however, in the case k = —1 (open universe), this inflexion point is located after 
the end of conformal time (r > Tf) and is therefore physically irrelevant. The existence 
of such a point is of interest only in the case of a spatially closed universe. Notice that in 
that case, the experimental observations tell us that the present value Tq of r (i.e., todays' 
date) is bigger than r/ (and, of course, smaller than Tf). 

The values Tf < t < Tg correspond to a classicaly forbidden region (negative dimen- 
sionless temperature). The values Tg < t < 2ujr correspond to universe in contraction, 
ending with a big crush . . . Finally, one should remember that T is a (doubly) periodic 
function, so that the same analysis can be performed in all the intervals [2puJr,2{p+ l)LOr]- 

Let us repeat that r is a conformal time and that t ^ oo when r — > Tf, so that, in 
both cases k = ±1 we are indeed in a ever expanding universe starting with a big bang. 
The curve obtained for k = has the same qualitative features, but for the fact that the 
inflection point moves to the end of conformal time: tj = Tf in that case. 

The shape of T(t), as given by figures I and I is in full agreement with the qualitative 
discussion given in the previous section. Such plots were already given in [|[ where Weier- 
strass functions had been numerically calculated from the algorithms obtained in the same 
reference. The same plots can now be obtained easily using for instance Mathematica, (the 
last versions of this program incorporate routines for the Weierstrass functions V, C and 
fj). Note, however, that intensive calculations involving such functions should probably 
make use of the extremely fast algorithms described in reference [11 1; these algorithms 



use, for the numerical evaluation of these functions, a duplication formula known for the 
Weierstrass V function. Moreover, evaluation of elliptic integrals (for instance caculation 
of periods in terms of invariants) can also be done by an iterative procedure (based on 
classical properties arithmetic-geometric mean), this is also described in reference pT[| . 
The duplication formula reads: 



4[4J'3(z)-9,P(z)-93] 

Therefore, one obtains a relation between the values of the reduced temperature T at 
(conformal) times r and 2t which, in the case k = +1, reads 

T(2t) — - = —2(T(t) — -I [(^(''") ~ 2)^ ~ j]^ 

^ ^ 2 ^^^2^ 4[4(r(r) - i)3 - 3(r(r) - i) - (1 - 2A)] 
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Figure 5: Evolution of the reduced temperature. Open case. 



When the universe is very young (r near zero), one may approximate the V{z) function 
by 1/z^ and therefore the reduced temperature T(r) by 

It may be useful to notice, at this point, that when the cosmological constant is zero 
(so that we have also A = for the reduced cosmological constant), one obtains 

r(r) 



2 sin^ f 



whereas, when the cosmological constant A is equal to the critical value Ac (so that A = 1), 
one obtains 

T(r) = 1 + ^ 

^ ^ 2sinh2 5 

As already discussed previously, recent measurements imply that A > 1, so that T(r) 



is given by formula 22 
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4.2.3 Not Neglecting Radiation : the Case a 7^ 

This case, technically a little bit more difficult, is detailed in Appendix 1. 

4.3 Determination of all Cosmological Parameters: a Fictitious Case 

Here we summarize the procedure that should (ideally!) be followed, in order to specify 
all cosmological parameters of interest. 

4.3.1 A Simple Procedure 

• From experimental results on the present values of the Hubble function Hq and of the 
temperature of the cosmic background radiation, calculate the radiation density 
Q° (use formula]^. 

• Improve the experimental bounds on the present values of matter density 17^^ and 
cosmological constant density J^a (so, without assuming that k is zero, of course). 

• From the "Cosmic triangle relation" (eq. obtain the present value of the curvature 
density function 17^. The sign of this quantity determines also the opposite of the 
constant k = ±1 (or 0). 

• From these results on densities f]'s, compute the dimensionless constant parameters 
a and A (use equations 13 and 12). 



From a and A compute the two Weierstrass invariants (72 and (73 that characterize 



the elliptic curve associated with our universe (use equations IC and O) 



The present value of To = T(ro) of the reduced temperature is obtained from equation 
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If one decides to neglect radiation (a ~ 0), one can plot directly the reduced 
temperature in terms of the conformal time T{t) = Q['P{t; §2, gs) + 1/12], where 
V is the Weierstrass elliptic function, by using, for instance, Mathematica. The 
present value Tq of conformal time is obtained by solving numerically the equation 
To = Q['P{to; g2, gs) + 1/12]. Other interesting values like rj (the end of confor- 
mal time) or tj (inflection point, only interesting if A; = +1) can be determined 
numerically; for instance rj is obtained, with Mathematica, thanks to the function 
FindRoot (remember that r/ is the first positive zero of T{t)). The value of the 
half-periof uJr can be found by evaluation of an elliptic integral but it is simpler 
to determine it by solving numerically the equation = since t = oJr corre- 

sponds to a global minimum. Another possibility, still with Mathematica, is to use 
the function W eierstr ass Half Per iods[{g2, gs}] (this function returns a set of two 
elementary periods, and twice our tOr coincides with the real part of one of them). 

If one decides to keep a 7^ 0, one has to use the results of Appendix 1 and determine 
Tf and Tg first; these two real zeros of T(r) (within a periodicity cell) are given by 
the elliptic integrals given in Appendix 1; they have to be evaluated numerically (the 
value Tg was just equal to 2uJr — tj in the case with no radiation, but here, it has 
to be determined separately). The dimensionless temperature T(r) is then given, in 
terms of the Weierstrass (" function, by formula 

In all cases, the present value of the deceleration function (a bad name since it is 



experimentally negative!) Qo is given by equation |15 
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Figure 6: Evolution of the reduced temperature in a neighborhood of To 



The value of Hubble constant allows one to determine all the dimensional quantities, 
in particular the critical length parameter Ac (formula |l^ , the cosmological constant 
A (= AAc), and the present values of the the cosmic scale (or "radius") Rq (formula 
llH) and of the cosmic time to (formula |l9|) . 



4.3.2 An Example 

Here we follow the previous procedure assuming fully determined values for the Hubble 
constant and density parameters and J7a- Of course, such precise values are not 
yet experimentally available and we have to make a random choice (compatible with 
observational bounds) in order to illustrate the previous "simple procedure" leading to 
the determination of all cosmological parameters of interest. In other words, what we are 
describing here is only a possible scenario. For simplicity reasons; we shall neglect the 
influence of radiation (so a ~ 0). 

We take experimental densities Vtm = 0.4 and JIa = 0.8. Then = —0.2 and 
k = +1. The reduced cosmological constant is A = 108. The Weierstrass invariants are 
52 = 1/12 and = —1.00463. The evolution of the dimensionless temperature is T{t) = 
6V{t, {(72, 93}) + 1/12 and its present value turns out to be equal to To = 3; this is a curious 
coincidence (see section 3.3) since To has no reason, a priori , to be equal to the temperature 
of the cosmological black body radiation Tg. The present value of the conformal time, 
obtained by solving numerically the equation T{to) = To is Tq = 1.369. The end of 
conformal time, obtained by solving numerically the equation T{Tf ) = is rj = 1.858. 
The (real) half-period cOr of the function T(r) may be obtained by solving numerically the 
equation T{t)' = is = 2.662. The conformal time for which T(r) has an inflexion point 
is obtained by solving numerically the equation T[t)" = and is r/ = 1.691. Notice that in 
that particular universe, we have To < tj (< Tf), so that we have not reached the inflexion 
point, yet. Using now the value of the Hubble constant Ho = (13.25 x 10^'')^^c?7i~^, 
we find a characteristic length scale of matter equal to Ac = 1.266 x 10^^^ cm^^ so that 
the cosmological constant itself is equal to A = 1.367 x 10~^^c?n^^. Finally we find 
the cosmic scale Ro = R{to) = 2.96 x 10^*^ cm = 3.14 x IQi^^yrs and the cosmic time 
to = t{To) = 1.24 X 10^*^ cm. Figure ^ gives the behaviour of T(t) in a neighborhood of 
To In this example we neglected the influence of radiation described by the dimensionless 
parameter a (given by equation [l^) ; using the experimental value of the density parameter 

= 5.1 X 10~^, one finds a = 2.77 x 10~^. This does not modify Tq but modifies the 
values of the Weierstrass invariants §2 and (73; also, the behaviour of T(r) is not the same. 
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Taking this value into account together with the results of Appendix 1 leads to a slight 
modification of the values of Tq and therefore of to- 
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5 INFLUENCE OF THE COSMOLOGICAL CONSTANT 
ON THE REDSHIFT. LARGE SCALE STRUCTURES 
AND GEOMETRY OF THE UNIVERSE. 

5.1 Cosmological Constant Dependence of the Redshift Function 

We shall now discuss the redshift, z, as function of To and of the difference 

5 = To-T 

where r is parameter time (conformal time) r of the emitter and parameter time, and To 
the parameter time of the observer. The redshift is given by 

, _ Rjro) _ T{t) 
R{t) T{to) 

Assuming, as before, that we are in the situation A > 1,Q! = 0, which is compatible with 



recent observations, we obtain immediately, from equation 22, the expression 



z = C[V{to -5) + ^] -I (23) 

with 

C = [V{to) + (24) 

Remember that the Weierstrass elliptic function V is characterized by the two Weierstrass 
invariants 92 = j2 ^t^'^ 93 = ^(^^ ~ '^^)- order to illustrate the influence of the 
cosmological constant on the behaviour of the redshift function, we continue our previous 
example (section |4.3.2| ), therefore taking fc = +1 (a spatially closed universe). To = 1.3694, 
A = 108 and plot z as a function of 5 (see figure 0). The value S = To (i.e., r = 0) 
corresponds to a photon that would have been emitted at the Big Bang whereas 6 = (i.e., 
r = To) corresponds to a photon emitted by the observer himself (no redshift). The dashed 
curve (the lower curve) gives the corresponding redshift when the cosmological constant 
is absent (A = 0). It is clear from these two curves that the influence of the cosmological 
constant becomes stronger and stronger when the geodesic distance 6 increases. 



5.2 Observer Dependence of Redshift Values and Large Scale Geometry 
of the universe 

As already announced in the introduction, we now want to consider the problem of com- 
parison of "redshift charts": Given two observers, far apart in the universe, how do we 
compare the redshift that they will record if they look at the same astronomical object ? 

Actually, one may consider (at least) two very different situations. The first is to 
assume that the two observers perform their measurement at the same time (a meaningful 
notion since we are in a homogeneous space-time), in that case, it is clear that the photons 
that they are receiving were not emitted at the same time by the cosmological source 
that they are studying. The other situation is to assume that both are interested in the 
same cosmological event (a star belonging to a distant galaxy turning to a supernova, for 
example) and our two observers record the redshifted light coming from this event when 
the light reaches them (but in general this light will not reach them at the same time). 
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Figure 7: Cosmological constant dependence of the redshift function 

5.2.1 Comparison of Measurements Made at the Same Time 

We first suppose that both observers perform their observation at the same time. It is 
natural to specify events by a pair (a, S) where S" is a point in space (Sun for instance) 
and a a particular value of the conformal time of the universe. Astronomers on earth 
(S), nowadays (tq), may decide to observe and record the redshift z of an astronomical 
object X (a quasar, say). At the same time Tq, astronomers located in the neighborhood 
of a star P belonging to a distant galaxy may decide to observe the same astronomical 
object X and record its redshift Z. The problem is to compare z and Z. Of course both 
z and Z are given by the previous general formula expressing redshift as a function of 
the time difference between emission and reception, but this value is not the same for S 
and for P. In order to proceed, we need a brief discussion involving the geometry of our 
three-dimensional space manifold. 

For definitness, we suppose that we are in the closed case (A; = 1), i.e., the spatial 
universe is a three-sphere S*'^ . The reader will have no difficulties to generalize our formulae 
to the open case (essentially by replacing trigonometric functions by the corresponding 
hyperbolic ones). The fact that universe expands is taken into account by the evolution 
of the reduced temperature as a function of conformal time r and, as far as geometry is 
concerned, wc may analyse the situation on a fixed three-sphere of radius 1. The choice of 
the point P - wc shall call it the "Pole" but it is an arbitrary point - allows one to define a 
notion of equator and of cosmic latitude: the equator, with respect to P, is the two-sphere 
(a usual sphere) of maximal radius, centered on P and the cosmic latitude £{S) of the Sun 
S is just the length of the arc of great circle between S and the equator; this great circle 
is the geodesic going through the two points P and S (that we suppose not antipodal!). 
The astronomical object X under simultaneous study of P and S is also characterized by 
a cosmic latitude i{X). It is obvious that £{P) = |. It is also clear that the (conformal) 
time difference S = To — t between observation (by S at time Tg) and emission of light 
(by X at time r) is nothing else than a measure of the arc of geodesic defined by the two 
points S and X on the unit three-sphere. We have therefore a "triangle" PXS whose 
"edges" are arcs of great circles and the lengths of these three arcs are & = f — ^(3) 
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(from P to S), a = ^ — £{X) (from P to X) and c = 6 (from S to X). The last piece 
of information that we need is a measure a of the angleQ, as seen from Sun (S) between 
the direction of P and the direction of X. It is a priori clear (draw a triangle !) that 
there is a relation between i{S), i{X), 5 and a, a relation that generalizes the well known 
formula = 6^ + — 26c cos a valid for an arbitrary triangle in euclidean space. Here is 
the formula that we need: 

sin ^(X) = cos 5 sin £{S) — sin 5 cos a cos £{S) (25) 

The proof of this formula is given in Appendix 2 and uses the fact the three-sphere 
carries a group structure: it can be identified with SU{2) or with the unit sphere in the 
non-commutative field of quaternions H. The uninterested reader may take the above 
formula for granted but the technique used in our proof is of independent interest and 
may be used by the reader to solve problems of similar nature. 

The formula ^ gives in particular the (conformal) time difference 6 between S and X, 
when all the points of our spatial universe are characterized in terms of cosmic latitude 
with respect to an arbitrary reference point P (the pole). The equation giving 5 should 
be solved numerically in general, but notice that the equation simplifies, and give rises to 
an analytic expression, when we choose X on the equatorial two-sphere defined by P. In 
that particular case, £{X) = and 

6 = arctan (26) 

cos(a) 

Coming back to our problem of comparing redshifts, we find that the redshift of the 
object X, as observed by P at conformal time Tq, is 

z = c[v{To-{7T/2-e{x))) + ^]-i 



with C still given by equation 24, whereas the redshift of the same object X, as observed 



by S (Sun) at the same time Tq is given by equation 53, with 5 determined by equation 



25| ; as before, £{X) and £{S) refer to the cosmic latitude of X and S with respect to the 
reference point P (pole) and a is the angle between the sighting directions SP and SX. 

In the particular case of an astronomical object X belonging to the two-sphere which 
is equatorial with respect to the reference point P, the formulae can be simplified: The 
redshift of X as, recorded by P is then 



_ [P(r, - vr/2) + ^] 



whereas, as recorded by S", it is 



[Viro - arctan(^^)) + 



1 (28) 
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The redshift of equatorial objects X, as seen from P, can be numerically computed from 
equation and it is a direction-independent quantity. The redshift of the same equatorial 
objects, as seen from S is clearly direction dependent. The existence of a direction- 
independent significative gap in the distribution of quasars, as seen from a particular point 

^This symbol a has of course nothing to do with the one that was previously used to denote the radiation 
parameter 



23 



P of the Universe, around the value Z of the redshift, would certainly be an example of 
a (remarkable) large scale structure, but such a gap would become direction dependent 
as seen from another point S (the Sun); the previously given formulae would then be 
necessary to perform the necessary change of redshift charts allowing one to recognize the 
existence of these features. 

Notice however that, in this particular example {X equatorial with respect to P) we 
have to suppose Tq > tt/2 since, in the opposite case, the light coming from X cannot be 
recorded by P. Such an hypothesis is not necessarily satisfied and actually is not satisfied 
when we use the values given in our numerical example of section 4.3.2 since we found, 
in that case, Tq = 1.369 < 7r/2. In such a situation, if we want to compare measurements 
between S and P made at the same time, we have to look at an object X which is not 
equatorial with respect to P but is such that its light can reach both S and P at time 
To] the first condition (observability from P) gives tt/2 — i{X) < Tq i.e., i{X) > ir/l — Tq 
whereas the second one (observability from S) reads S < To, but 6 is still given by ^ so 
that this condition translates as a condition on the angle a between the sighting directions 
of X and of P as seen from S : this angle should be small enough with a maximal value 



obtained by replacing 6 by Tq in equation 25 and solving for a. 



Another possibility to illustrate the above formulae is to compare the redshifts of X 
measured at time Tq from two positions Si and 5*2 that should have small enough cosmic 
latitudes since the two conditions on 5i and 82 (the geodesic distances between X and 
these two positions) are 61 < Tq and 62 < Tq. To simplify the calculation we suppose that 
Si and S2 belong to the same meridian (going through P) of our spatial hypersphere, 
and call ai (resp. 02) the angle made between the sighting directions of X and of P, 
as seen from Si (resp. from S2)- Using the technique explained in Appendix 2 (write 
S2X = S2S1S1X), the reader will have no difficulty to prove the following formula that 
generalizes equation |2^: 

62 = arccos(cos(^(S'2) —£{Si))cos6i — sin (^(5*2) — ^(^i)) sin cos ai (29) 

To simplify even further the calculation, we take X equatorial with respect to the reference 



point P; therefore the previous inequalities on 61 and 62, using equation 12^, give, in turn, 
the following two conditions on the cosmic latitudes of the two observers: tan ^(5*1) < 
tan To cos «! and tan^(52) < tanroCosa2- The redshift zi of X as measured by 5*1 is 



still given by equation E3^ but we can use the particular equation 26, valid for equatorial 



objects, therefore zi = C\P{to — 5i) + j^] — 1 with 61 = arctan(^^^^^^) whereas Z2 
C[V{to — 62) + ^] — 1 and 62 given by equation 29. 



In order to illustrate these results in the study of large scale structures, we choose, 
for the cosmological parameters, the particular values given in section 4.3.2; furthermore 



we fix the values of ^(5*1) (take tt/6) and of ai (take 7r/4). Then 5i, is given by 26, and 
the redshift of X, as measured by 5*1 is zi = 3.42. The redshift Z2 of X, as measured 
by 5*2 is given by the above formula, but since everything else is fixed, it becomes a 
function of the cosmic latitude ^(6*2) only. For definitness we choose P, Si and 5*2 in 
the same hemisphere. The following figure (P) gives Z2 as a function of ^(^2) in the range 
[tt/6 — 7r/12, 7r/6 + 7r/12]; of course this curve intersects the horizontal line zi = 3.42 when 
£{S2)=£{Si) = tt/6. 



^existence of such a gap was investigated in references 
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Figure 8: Dependence of the redshift function on the cosmic latitude of an observer 



5.2.2 Comparison of Measurements Relative to a Single Event 

We now examine the situation where both observers (take P and S) are interested in the 
same cosmological event (a star belonging to a distant galaxy turning to a supernova, for 
example) and record the redshifted light coming from this event when the light reaches 
them (both measurements are not usually performed at the same time). Here we do not 
have to assume that Tq > vr/2, even if X is equatorial with respect to P since the light 
emitted by X will reach P anyway, in some future. Keeping the same notations as before, 
the cosmological event of interest took place at conformal time Tq — 6, and the light will 
reach P at time Tq — 6 + it/2 — £{X). The redshift measured by S is, as before, given by 



equation 23 



but the redshift measured by P is 



[Piro) + ^ 



[V{To-6 + 7T/2-iiX)) + i-,] 
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A The Case With Radiation {a ^ 0). Analytic Solution. 



General features When a ^ 0, the two poles of T(r) are distinct (and each of them 
is of first order). The three roots ei, 62, 63 of the equation P{y) = (^l|) add up to zero. 
They are given by ± i^iA - B) and A + B, where 

A = l/2[g3 + v /^a727 ]^/=^ (30) 
B = l/2[53 - \/^a727]1/3 (31) 

and where A is the discriminant 

A = 5i - 27 gl = 3-3a3A(A - A+)(A - A_) 

with 

A± = + Uka + 1 ± (8ka + if/^) 

32a'^ 

As a — > 0, we have 

A^2-^3-3A(1-A) 
indeed, A+ — > 00 and A_ — > 1 and, to first order in a, we have 

A, Y^(1 + 12M 

A- ~ k{l - 3a/fc) 

If A = - something that is nowadays ruled-out experimentally but that used to be called 
the "physical case" ! - two zeroes of P{t) coincide (then also, two zeroes of Q{t)); in 
this situation one of the two periods becomes infinite and the elliptic function degenerates 
to trigonometric or hyperbolic functions. The analytic study, in that situation, is well 
known. 

We know assume that a 7^ and A_ < A < A-|- This is the situation that seems to 
be in agreement with the recent experiments (see the discussion given previously). For 
such values of A, the polynomial Q{T) has two negative roots Ta < Th < (see fig ||). 
Here 53 < and the discriminant A is (strictly) negative, as well. The three roots of the 
equation P{y) = can be written 

ei = a — ib 
62 = —2a 
es = a + ib 

with a > and 6 > 0. 

In the case a 7^ 0, it is of course still possible to express T in terms of y (using the 
fractionally linear transformation with y = 'P{t; g^, gs), the Weierstrass "P-function. 
This is actually not very convenient and it is much better to use either the Weiertrass 
("-function (^(r) or the Weierstrass fi-function (t(t). 

The function (^(r) - which is an odd function of r is not elliptic but it is defined by 

and the requirement that 1/t — ((r) vanishes at r = 0. Although is not elliptic, the 
difference C{t — u) — C(t — v) is elliptic of order 2, with poles at u and v, when u ^ v 
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Figure 9: Typical evolution of the reduced temperature (case with radiation) 



are two arbitrary complex numbers. In terms of the quantities R{t) and T(r) read 
immediately: 

^Rir) = = - + C(r,) - C(r - y) - C(r;) (32) 

where the two zeros of T{t), < tj < Tg < 2^,. are given by (use figure |l|) : 

dT 

V 



Tf + 2 



dT 



In VW) 

The quantities and Th are defined on the graph given on figure they are numerically 
determined by solving the equation Vak{T) = ^ (for instance use the function FindRoot 
of Mathematica) . The curve T{t), in the closed case, has the shape given by figure ^ in 
full agreement with the qualitative discussion given in section 4; of course, only the first 
branch of this curve, from the big bang (r = 0) to the first zero Tf of T{t) (the end of 
conformal time) is physically relevant. It is instructive to compare this function with with 
the a = case (no radiation), given in figure ^. In both cases, we see that the function 
r(r) is negative between tj and Tg (its next zero), but here, it becomes infinite again for 
a value rj + Tg of conformal time which is strictly smaller than 2L0r ; after this singularity, 
we have a negative branch (a novel feature of the case with radiation) in the interval 
]r/ + Tg,2wr[; the period, on the real axis is, as usual, denoted by 2u}r- One shows that 
2ojr = Tf ~^ Tg + Tc where 

^'^ dT 

2 



Such plots were already given in Q where Weierstrass functions had been numerically 
calculated from the algorithms obtained in the same reference. 

The same reduced temperature can be expressed in terms of the Weierstrass function 
ct(t). The even function ct(t) is not elliptic either and is defined by 

I da 
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and the requirement that (T(r) should be an entire function vanishing at r = 0. Although 
not elliptic, the quotient is elliptic of order 2, with poles at U3, U4 and zeroes 

at ui, U2, whenever ui, U2, W3, 1x4 are complex numbers such that ui + U2 = + U4 . In 
terms of a, the quantity T{t) reads immediately: 

a{T)a{T -Tf- Tg) 

where the constant c is given by 
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B Geometrical Relations in a Geodesic Triangle 



The unit three-sphere will be here identified with the unit sphere in the non-commutative 
field of quaternions M ~ M^. To each X e C W^, X = {Xq, Xi, X2, X3) we associate 
the quaternion 

X = Xo+Xii + X2j + X3k (33) 

We shall use the non commutative multiplication rules = = A;^ = — 1, ij = —ji = k, 
jk = —kj = i and ki = —ik = j. The norm square of X is given by |Xp = XX = Xq + 
Xf + X2 + X^ where X is the quaternionic conjugate of X ( with i = —i, j = —j, k = —k). 

It is convenient to write X in a form analoguous to the representation X = e^^ = 
cos 5 + i sin (5 of a complex number of unit norm; in this elementary situation, the angle 
5 is the geodesic distance (arc length) between the unit 1 and the complex X, whereas 
i (unit norm) can be considered as the tangent (sighting) direction from 1 to X along 
an arc of geodesic (since i represents the vector (0, 1) tangent to the unit circle at the 
point (1,0)). In the three dimensional situation, the unit 1 is a particular point of the 
three-sphere 5^ and we shall write, in the same way 

X = e^^ = cos5 + Xsin5 (34) 

where 6 is the geodesic distance (arc length) between the unit 1 and the quaternion X, 
while X is a quaternion of square —1 and unit norm representing the three dimensional 
vector tangent to at the point 1 in the direction Y (sighting direction). More precisely, 
let X^, £ = 1,2, 3, be this three dimensional unit vector (Xf + X| -|- X| = 1) and call 

X = Xii + X2k + X^k 



Then XX = — 1 and XX = -|-1. Comparing equations 33 and 34 gives 



1. = ^ (35) 
smo 

and 

cos 5 = Xq (36) 

Take Tt and if two three-dimensional vectors (not necessarily of unit norm), call 
u = uii + U2j + u^k, V = vii + V2j + v^k the corresponding quaternions and u, v the 
corresponding normalized unit vectors. An elementary calculation, using the previous 
multiplicative rules, leads to uv = — cosq -|- 1? x if, but \lf x if] = \lf\ |lf | sina, so 
that 

uv = \lf\ \ lf\{ — cos a + if X if sina) 

where a is the angle between vectors lf,lf, and I?' x if is the normalized vector product 
of these two three-dimensional vectors. For normalized vectors, we have in particular 
{uv)o = — COSQ, or better, 

{uv)o = +cosa (37) 

Now, let P (a reference point "Pole"), S (Sun) and X (a quasar) three points of 
S''^ C HI ~ that we represent by quaternions also denoted by P, S and Y. We call i{S), 
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resp. i{X), the latitude of S (resp. of X) with respect to P, counted positively whenever 
P and S (resp. X) are in the same hemisphere, and counted negatively otherwise. 

Let Y = SX. From equations ^ and 36, we find that the geodesic distance 5 between 
X and S (or between 1 and Y) is such that cos 5 = Yq so, 

6 = arccos{SX)o (38) 

and that the sighting direction from S" to X is given by a unit vector Y with three 
components 

Y, = ^ (39) 
sm 

We have _ 

Y = SX = cos 6 + Y sin 6 = cos 6 + SX sin 6 

Notice that, in the same way, 

PS = cos(^ - i{S)) + sin(| - e{S))TS = sm{i{S)) + cos{£{S))TS 

and 

PX = sm{£{X)) + cos{i{X))PX 

Using the fact that SS = 1, we remark that 

PX = PSSX 



This last equality provides the clue relating the three sides of our geodesic triangle. It 
implies in particular (PX)o = {PSSX)q. The left hand side is given by (PX)q = sm{i{X)) 
and the right hand side comes from 

(PS)(SX) = {sm{£{S)) + cos{£{S))TS){cos6 + sm6fx) 
= sin £(S) cos 5 + cos £{S) sin 6PS SX + 
sme{S) sm6SX + cosi{S) cos6PS 
{(PS)(SX))o = sm£{S) cos6 + cos£{S) sm6(Psfx)o 

By equation ^ we know that the cosine of the angle a between vectors SP and Sjl, 
respectively represented by the quaternions SP and SX, is given by the 0-component of 

the quaternionic product {SP){SX) = {PS){SX). 

To conclude, one finds the following formula relating the lengths of the three sides of 
the geodesic triangle SXP, together with the angle a: 

sm{£{X)) = sin^(S') cos 5 — cos£{S) sin 5 cos a 
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